12.1 FortranとCの混合
Fortran ソースコードとCのコードを混在させる場合でFortran 2003のC相互運用を利用しない場合には、 以下の点に留意してください。
12.1.1 ネーミング規約
外部手続き名と共通ブロック名は小文字にアンダスコアを付加したものです(これは 標準UNIX f77と同じです)。メインプログラム単位の名前は‘main’と
なります。
モジュール変数とモジュール手続き名はモジュール名(小文字)に
‘_MP_’を付加、さらに言語要素名(小文字)を付加する形で構成され
ます。
注:これらの規約は-compatible オプションを使用した場合には異なる可能性があります。 Windowsにおいてこのオプションを指定すると、 外部手続き名と共通ブロック名は大文字でかつアンダスコアが付かないものとなる一方、 モジュール要素はモジュール名を小文字ではなく大文字で持つことになります。 (64ビットシステムでは常に-compatibleです。) その他の多くの環境では、 このオプションの使用により名前にアンダースコアが既に含まれている場合に、 (通常付け加えられるもの後に更に)アンダースコアが追加されます。
12.1.2 初期化と終了
メインプログラムがFortranで書かれていない場合には初期化ルーチン
__nAGf90_rts_init(int argc,char*argv[])
または入出力初期化ルーチン
__nAGf90_io_init(void)
のいずれかをCからコールする必要があります。
__nAGf90_rts_initルーチンはFortranの浮動小数点環境を初期化すると共に、
Fortran 2003の組込関数及びF90_UNIXモジュールからのコマンド行引数に対するアクセスを可能にします。
またプログラムは__nAGf90_finish(int status)によって終了させる必要がありま
すが、通常のC終了ルーチンを用いてすべてのFortranファイルのクローズ、Fortran
入出力バッファのフラッシュを行う前に、__nAGf90_io_finish(void)をコールする
というアプローチも取れます。(フラッシュを行わないとオープンされている出力ファイルが壊れる可能性があります)
12.1.3 呼出し規約
以下のセクションでは、nAG Fortranコンパイラを用いてコンパイルされたFortran プログラムによって用いられる呼出し規約を、Cの用語を用いて詳細に記述します。 従ってこの情報はFortranとCのコードを混在させたいと考えている人にとって最も 有用と言えます。
近代FortranのFortran 77サブセットを対象とした規約はデファクトである
Unix f77規約と互換性があります(-f77 オプションを用いずに
コンパイルされたCOMPLEX関数を除く)。
12.1.4 データタイプ
CからFortranへの交信に役立つデータタイプの定義は、通常/usr/local/lib/nAG_Fortran(UNIX系のシステムの場合)に置かれているdope.hとnagfortran.hファイルの中にあります。
| Fortranデータタイプ | Fortran精度 | C typedef名 |
| INTEGER | 8 bits | Integer1 |
| 16 bits | Integer2 | |
| 32 bits | Integer または Integer3 | |
| 64 bits | Integer4 | |
| LOGICAL | 8 bits | Logical1 |
| 16 bits | Logical2 | |
| 32 bits | Logical | |
| 64 bits | Logical4 | |
| REAL | single | Real |
| double | Double | |
| double-double | DDReal | |
| quadruple | QReal | |
| COMPLEX | single | Complex |
| double | DComplex | |
| double-double | DDComplex | |
| quadruple | QComplex |
DDRealとQRealはほとんどのシステムで同一です。ただしSun Solarisにおいてはこれら は異なります(後者はIEEE互換の128ビット浮動小数点型となります)。
12.1.4.1 ポインタ
スカラで多相でなく文字型でもないPOINTER型は単にオブジェクトに対するCポインタです。
配列POINTER型は対象の配列を記述するドープベクトル(dope vectors)です。
単なるアドレスとは異なり、これらのドープベクトルは連続な配列中における不連続
な部分配列を直接記述することができます。
詳細については下記(DopeN とChDopeN)をご参照ください。
多相のドープベクトルはCSDopeNであるCLASS(*)ポインタを除き NPDopeN です。
12.1.4.2 構造型
Fortranの構造型はCのstructに翻訳されます。
BIND(C)またはSEQUENCEが使用された場合には、構造体内の項目の順序は
構造型定義中における順序と同一になります。
そうでない場合には、並び方を改善し記憶容量を減らすべく順序の入れ替えが行われ、
より大きな型が構造体の先頭の方に置かれることになります。順序の確認はCの出力コードを調べる事で行えます。
12.1.4.3 サポートタイプ
Charunsigned char
基本(単一バイト)文字記憶のデータ型Char2unsigned short
16ビット(JIS及びUCS-2)文字記憶のデータ型Char4unsigned int
32ビット(UCS-4)文字記憶のデータ型Substringstruct { Char *addr; Chrlen len;}
単一バイト(基本)CHARACTER列を記述します。 無指定文字長の基本CHARACTER変数に使用され、 可変長のスカラnon-POINTERCHARACTER関数、 及びすべてのPOINTER CHARACTER関数の戻りタイプとして使用されます。Substring2struct { Char2 *addr; Chrlen len;}
16ビット(JISもしくはUCS-2)CHARACTER文字列を説明し、 無指定文字長の16ビットCHARACTER変数、 及び可変長スカラ非POINTERの16ビットCHARACTER関数の戻り型、 及び全てのPOINTER16ビットCHARACTER関数で利用されます。Substring4struct { Char4 *addr; Chrlen len;}
32ビット ISO_10646(USC-4)CHARACTER文字列を説明し、 無指定文字長の32ビットCHARACTER変数、 及び可変長スカラ非POINTERのUSC-4CHARACTER関数の戻り型、 及び全てのPOINTERUSC-4CHARACTER関数で利用されます。Offsetint,longまたはlong long
アドレス付けと添え字計算のための整数型です。これは32ビットシステム及び小型64ビット システム上ではint(32ビット)、大型の64ビットシステム上では64ビット整数 型です。Chrlen- 通常
int(64ビットのWindowsではlong long)
文字長を表すための整数型です。 Pointer- 通常
char *
任意の型を参照するために、及びポインタ算術用に使用されるバイトポインタです。 Tripletstruct { Offset extent,mult,lower;}
配列次元のパラメータを含みます。extentはその次元のサイズを、multは 刻み幅(連続した要素間の距離(単位バイト))を、lowerは下限を表します。 これはDopeN,ChDopeN構造体の成分です。DopeNstruct { Pointer addr; Offset offset; Triplet dim[N];}
すべての非多相non-CHARACTER配列(構造型の配列を含む)に対するドープベク トルです。N は次元数を意味し1から7までです。addrは最初の要素のアド レスであり、dimはそれぞれの次元を記述します。offsetは添え字計算 に足されるオフセット、すなわちSUM(mult*lower)を表します。これはPOINTER配列関数に対する戻りタイプとして使用されます。それに対するポイン タは形状引継ぎ配列、POINTER配列引数用の引数タイプとして使用されます。ヌル化された配列ポインタはヌルポインタである
addrフィールドを持ちます。 ゼロサイズの配列はヌルポインタではないaddrフィールドを持つので注意し てください。ChDopeNstruct { Pointer addr; Chrlen len; Offset offset; Triplet dim[N];}
DopeN 構造体にCHARACTER長を指定するlen成分が加わった ものです。ArrayTemp_typestruct { type *addr; Offset extent[7];}
type がInteger1,Integer2,Integer,Integer4,Logical1,Logical2,Logical,Logical4,Real,Double,QReal,Complex,DComplex,QComplexのいずれかである連続配列を記述します。 それはnon-POINTER配列関数に対する戻りタイプとして使用されます。 配列の次元数より後のextent値は不定である点に注意してください。ArrayTemp_Characterstruct { Char *addr; Chrlen len; Offset extent[7];}
連続なCHARACTER配列を記述します。他の配列型と同じですが、CHARACTER長に対するlen成分が追加されています。ArrayTemp_DerivedArrayTemp_Characterと同義
任意の構造型の連続配列を記述します。 この場合のlenは構造型配列要素のサイズを表します。
12.1.5 SUBROUTINE戻りタイプ
12.1.5.1 ラベル引数を持つSUBROUTINE
戻りタイプはintです。その値は制御を移す先のラベルに対するインデックス
(最初のラベルであれば1、以下同様)を表します。
ゼロ、または範囲外の値は制御の移行を伴わないことを示します。
12.1.5.2 ラベル引数を持たないSUBROUTINE
戻りタイプはvoidです。
12.1.6 FUNCTION戻りタイプ
12.1.6.1
スカラ、non-POINTER関数
-
INTEGER,LOGICAL,REAL
上記組込み型。 -
COMPLEX, -compatible オプション使用せず(もしくはWindows 64ビットモードではない)
精度に応じてComplex、またはDComplex。 -
COMPLEX, -compatible オプション使用(もしくはWindows 64ビットモード)
戻りタイプはvoidです。 型がComplexまたはDComplexの一時データのアドレスが関数の初期引数 として渡されます(結果はその位置に書き込まれます)。 -
CHARACTER, 固定長または長さ引継ぎ
戻りタイプはvoidです。 2つの追加の初期引数が関数に引き渡されます。1番目のもの(Char*,Char2*,Char4*のいずれか)は結果 を書き込む場所を示すアドレスを、2番目のもの(Chrlen)は結果の長さ(呼 ばれた関数が長さ引継ぎ型の場合)を示します。 -
CHARACTER, 可変長
戻りタイプはSubstring,Substring2,Substring4(上記参照)のいずれかです。 呼ばれた関数は文字列用の記憶領域を割付けます。その領域が必要なくなったときに それを解放するのは呼出し元の責任です。 -
構造型
構造体。
12.1.6.2
スカラPOINTER関数
すべてのPOINTER値関数に関し、結果がポイントする記憶領域は呼ばれる関数
内で割付けられていることもあれば、結果が既存の記憶領域をポイントしていること
もある点に注意してください。
-
INTEGER,LOGICAL,REAL,COMPLEX
適切な組込み型(例えばComplex*)へのポインタ。 -
CHARACTER
戻りタイプはSubstring,Substring2,Substring4のいずれかです。 -
構造型
構造体へのポインタ。
12.1.6.3
配列non-POINTER関数
-
組込み型
上記組込み型に対する適切なArrayTemp_構造体。 -
構造型
そのlen成分が構造体のサイズにセットされた形でArrayTemp_Derivedが 応答として返されます。
12.1.6.4
配列POINTER関数
-
CLASS(*)
配列の次元によりCSDope1,CSDope2, ...CSDope7のいずれか -
CLASS(derived-type-name)配列の次元によりNPDope1,NPDope2, ...NPDope7のいずれか -
非多相の非
CHARACTER型
配列の次元数に応じてDope1,Dope2, ... ,Dope7のいずれか。 -
CHARACTER
配列の次元数に応じてChDope1,ChDope2, ... ,ChDope7のい ずれか。
非多相の構造型配列はDopeN 構造体として返される点に注意してください。
12.1.7 引数タイプ
12.1.7.1
CHARACTER型
CHARACTER型の通常の引数はすべて、CHARACTER、マルチバイトCHARACTERを問わず、
CHARACTER言語要素の長さを示すChrlen引数を伴っています。この追加引数はすべての通常の引数よりも後、
引数並びの末尾で引き渡されます。CHARACTER引数が複数ある場合、その長さは左から右に順に引き渡されます。
これについて、32ビットのWindowsの場合で-compatibleオプションを指定した場合を除き、追加の引数はCHARACTERの直後に続きます。
このルールに対するもう一つの例外は形状引継ぎCHARACTER配列の場合ですが、この場合、仮引数の長さはドープベクトル中のフィールドから直接取得され、別個に引き渡されることはありません。
12.1.7.2
non-POINTER非ALLOCATABLEスカラ
-
non-
CHARACTER型
引数のアドレスが引き渡されます(例:INTEGER引数の場合のInteger*)。 -
CHARACTER
引数のアドレスが引き渡されると共に、引数の長さが別のChrlen引数として引数並びの末尾で引き渡されます。
12.1.7.3
POINTER及びALLOCATABLEスカラ
これらはポインタのアドレス(もしくは割付記述子)が引き渡される点を除くと、通常の場合と同様の形で引き渡されます(例:INTEGER POINTERの場合のInteger**)。
12.1.7.4
non-POINTER配列
-
引継ぎ形状
配列を記述する適切なドープベクトルのアドレスが引き渡されます (すなわち多相及び仮引数の型により、CSDopeN*,NPDopeN*,DopeN*,ChDopeN*のいずれか)。 配列は部分配列で記述できる限り連続である必要はありません。 -
その他
配列の最初の要素のアドレスが引き渡されます。CHARACTER配列の場合、それぞれの配列要素の長さはスカラに対する場合と 同様、引数並びの末尾でChrlenとして引き渡されます。 配列は連続でなくてはなりません。
12.1.7.5
POINTER配列
引継ぎ形状配列の場合と同様、適切なドープベクトルのアドレスが引き渡されます。
CHARACTER配列の場合、長さは別個の引数として末尾で引き渡されます。
12.1.7.6 手続き
手続きのアドレスが引き渡されます。CHARACTER関数の場合、長さは別個の引数として並びの末尾で引き渡されます。
関数が可変長型の場合、この長さは負の値となります。
12.1.7.7
OPTIONAL引数
OPTIONALな引数が.NOT.PRESENT()である場合、適切な型のヌルポイン
タが引き渡されます(例えばINTEGERスカラの場合、(Integer*)0が
引き渡されます)。
